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Abstract 

We present results of local spin density approximation (LSDA) pseudopoten- 
tial calculations for the perovskite structure oxide, bismuth manganite (BiMnOs). 
The origin of the differences between bismuth manganite and other perovskite 
manganites is determined by first calculating total energies and band structures 
of the high symmetry cubic phase, then sequentially lowering the magnetic and 
structural symmetry. Our results indicate that covalent bonding between bis- 
muth cations and oxygen anions stabilizes different magnetic and structural 
phases compared with the rare earth manganites. This is consistent with re- 
cent experimental results showing enhancement of charge ordering in doped 
bismuth manganite. 

1 Introduction 

Perovskite structure oxides exhibit a wide range of low temperature structural distor- 
tions associated with lattice instabilities of the prototype cubic structure (Figure |l|), 
including ferroelectric (e.g. PbTiC>3 and BaTiOa), antiferroelectric (e.g. PbZrOa) 
and antiferrodistortive (e.g. SrTi0 3 ). Both first- and second-order transitions are 
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observed, with a full spectrum of transition behavior ranging from displacive to order- 
disorder. The strong coupling between the electric polarization and the structural dis- 
tortions has led to the widespread use of perovskites such as lead zirconium titanate 
(PZT) in piezoelectric transducer and actuator applications]!]. 

In magnetic perovskites, different microscopic orientations of the spin-polarized 
ions give rise to different macroscopic magnetic symmetries. This results in rich 
phase diagrams, in which both the magnetic order and the structure depend on the 
temperature, pressure and chemical composition. A large amount of recent research 
activity has been focused on the rare earth perovskite-structure manganites, fol- 
lowing the observation of colossal magnetoresistance (CMR) in Ca-doped LaMnO3 0. 
During these recent studies, many rare earth perovskite manganites have been found 
to show strong coupling between their magnetic and structural, or magnetic and 
electronic order parameters. For example, a magnetically induced structural phase 
transition has been observed in Lao.83Sr .i7Mn03@] indicating strong coupling be- 
tween the local magnetic spin moments and the lattice structure. In Ndo.sSro.sMnOa, 
strong coupling between the magnetic spin moments and the electronic charge carriers 
was demonstrated when an electronic metal-insulator transition was induced by an 
external magnetic field ||. The large change in conductivity with applied magnetic 
field, which gives these materials potential technological importance as the read ele- 
ment in magnetic recording heads 0, is believed to originate from a similar type of 
phase transition |0]. 

Materials which have strong coupling between all three of the magnetic, electric 
and structural order parameters, resulting in simultaneous ferromagnetism, ferroelec- 
tricity and ferroelasticity, are known as multiferroics. Most multiferroic materials 
are complex structures with many atoms per formula unit, and more than one for- 
mula unit per unit cell. The large number of inter- ionic interactions has prevented 
isolation of the essential factors causing multiferroicity. As a result, the origin of 
multiferroism and the nature of the coupling between the magnetic, electric polar- 
ization, and structural order parameters are not well understood. In spite of the 
lack of fundamental understanding, a number of device applications have been sug- 
gested for multiferroic materials including multiple state memory elements, electric 
field controlled ferromagnetic resonance devices, and variable transducers with either 
magnetically modulated piezoelectricity or electrically modulated piezomagnetism|9|]. 
Also, the ability to couple with either the magnetic or the electric polarization offers 
an extra degree of freedom in the design of conventional actuators, transducers and 
storage devices. 

Bismuth manganite can be regarded as the "hydrogen atom" of multiferroics. 
Although information about BiMnC>3 is sparse||, indications are that it is simulta- 
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neously ferromagnetic and ferroelectric at low temperatures, and, because it has a 
simple structure, it is amenable to detailed study using first-principles techniques. 
Investigation of BiMnOs permits the study of the essentials of multiferroism with- 
out the problems associated with simulating compounds containing many atoms per 
formula unit. In addition, BiMnOs is strikingly different, both magnetically and 
structurally, from the rare earth perovskite manganites. BiMnOs is ferromagnetic 
with a triclinic structural distortion in its ground state, whereas the rare earth man- 
ganites are antiferromagnetic and either orthorhombic or hexagonal. The differences 
are unexpected because bismuth and the rare earths all form trivalent cations with 
similar ionic radii. The anomalous behavior of BiMnOs compared with conventional 
perovskite manganites could give rise to unusual and useful transport properties. 

The goal of the study described in this paper is to determine the origin of the dif- 
ferences between BiMnOs and the other perovskite manganites, both to understand 
the fundamental physics, and to assist in the optimization of perovskite manganite 
materials for novel device applications. To achieve this goal we evaluate electronic 
and magnetic properties of perovskite manganites using a plane wave pseudopoten- 
tial (PWPP) implementation of density functional theory (DFT) within the local 
spin density approximation (LSDA). Our choice of theoretical approach is influenced 



by successful first-principles studies of non-magnetic perovskite ferroelectrics[10. I I 



The origin of ferroelectricity in the prototypical ferroelectric perovskites PbTiOs and 
BaTi0 3 , and the reasons for the differences between them, were deduced using first 
principles density functional theory (DFT) calculations J12], [IT[] . It was shown in both 
cases that Ti 3d - O 2p hybridization is essential for ferroelectricity. In BaTiOs, the 
Ba 2+ ions are almost entirely ionically bonded, resulting in a rhombohedral struc- 
ture. In PbTiOs, however, the Pb 6s electrons are partially covalently bonded with 
the oxygen 2p electrons, favoring a tetragonal structure. 

Throughout this work, we use LaMnOs as our example of a typical rare earth 
manganite, and compare our calculated results for BiMnOs with those for LaMnOs. 
Since Bi is adjacent to Pb, near the right hand side of the periodic table, and La 
is adjacent to Ba on the left hand side, we anticipate by analogy with the earlier 
PbTi0 3 -BaTi0 3 study, that the high lying occupied Bi 6s electrons will contribute to 
the structural and magnetic differences between BiMn0 3 and LaMn0 3 . In addition 
we will find significant hybridization between the Bi 6p and O 2p electrons. 

The remainder of this paper is organized as follows: In Section we describe 
the plane wave pseudopotential LSDA implementation of density functional theory 
which we use for this study. Our PWPP implementation is shown to give results 
which are in excellent agreement with previously published all electron calculations 
for magnetic perovskite oxides. In Section [5] our results for BiMn0 3 and LaMn0 3 are 
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presented. Our results are summarized in Section f|. 



2 Computational Technique 

The calculations described in this work were performed using a plane wave pseu- 
dopotential implementation [13] of density functional theory fll4j] within the local spin 
density approximation. Plane wave basis sets offer many advantages in total energy 
calculations for solids, including completeness, an unbiased representation, and arbi- 
trarily good convergence accuracy. They also allow for straightforward mathematical 
formulation and implementation, which is essential in the calculation of Hellmann- 
Feynman forces [|15|] and in density functional perturbation theory calculations [|17 . 



However plane wave basis sets necessitate the use of pseudopotentials to model 
the electron-ion interaction, in order to avoid rapid oscillations of the valence wave- 
functions in the region around the ion cores. The accuracy and efficiency of ab 
initio pseudopotential calculations (compared with all electron calculations) is now 
well established for non-magnetic systems |T8[1 , but application of the method to spin- 
polarized magnetic systems is still in its infancy. Traditionally, magnetic materials 
containing transition metals have been studied using all electron methods with mixed 
basis sets, such as the linear augmented plane wave (LAPW)|19[, linear muffin tin 
orbital [2(| or Kor ringa- Kohn- Rost oker pTf approaches. 



Resistance to the use of the PWPP method to study these systems has been based 
on two perceived difficulties: 

1) First, the magnetic properties of the transition metals and their compounds 
originate from tightly bound ci-electrons, which might be expected to require a pro- 
hibitively large number of plane waves to expand their pseudopotentials. This prob- 
lem can be overcome by using the optimized pseudopotentials developed by Rappe et 
al. ||22|| . Optimized pseudopotentials minimize the kinetic energy in the high Fourier 
components of the pseudo wavefunction, leading to a corresponding reduction in the 
contribution of high Fourier components in the solid. The use of optimized pseudopo- 
tentials reduces the energy cut-off of the plane wave expansion in transition metal 
solids to around 60 Ry or less. 

2) Second, implicit within the pseudopotential approximation is the assumption 
that the exchange-correlation potential, V xc , is separable into a valence part and a 
core part: 

V xc (p, a) = [V xc (p, a) - V xc (p v , a v )] + V xc (p v , a v ). (1) 

Here p is the electron density, a is the spin polarization, and the superscript v refers 
to the contribution from the valence electrons. If the core electrons are spatially sep- 
arated from the valence electrons then this separation is rigorously correct. However 
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in transition metals, the (i-electrons have large wavefunction amplitudes near the core 
region. Therefore, because V xc is a non-linear function of both charge and spin polar- 
ization, use of Eqn. [TJ is inaccurate. Early pseudopotential calculations for magnetic 



materials attempted to circumvent this problem by using a small core radius |23], £4 



This resulted in a pseudo-wavefunction which was similar to the all electron wave- 
function and therefore could not be expanded efficiently using plane waves. A better 



solution is offered by the partial non-linear core correction scheme of Louie et al.[25 



in which the core charge density is explicitly retained in both the construction of the 
ionic pseudopotential, and in the DFT calculation. 

The PWPP method with optimized pseudopotentials and partial non-linear core 
corrections was first tested by Sasaki et al. for ferromagnetic nickel and ironpGf. 
The structural and magnetic properties calculated using the PWPP method were 
found to be in excellent agreement with results from all electron calculations. A 
similar study using different soft pseudopotentials obtained comparable results for the 
ground-state properties of nickel|27[]. More recently, Lewis and coworkers presented 
the first PWPP calculation for a magnetic compound (Cr02)f28|, and showed that the 
results agreed well with those of all electron LSDA calculations. Before presenting our 
results for BiMnC>3 in Section |3], we will compare the results of our PWPP calculations 
for a magnetic perovskite with published all-electron calculations |29| and show that 
similarly good agreement is obtained. 



2.1 Pseudopotential Construction 

For La and Bi, we constructed optimized scalar-relativistic pseudopotentials, using 
(for La) a 5s 2 5p e 5d l reference configuration with core radius (r c ) = 2.6 a.u., and 
(for Bi) a 6s°' 5 6p 1,2 6(i ' 3 reference configuration with r c = 3.0 a.u. In the La reference 
configuration the 4/ orbitals were unoccupied, and for Bi the filled 4/ shell was placed 
in the core. The pseudopotentials were optimized using 4 basis functions with a cutoff 
wavevector, q c , of 7.5 a.u. q c determines the convergence of the kinetic energy with 
respect to the plane wave cutoff energy in reciprocal space calculations. The resulting 
optimized pseudopotentials are shown in Figures ^| and The pseudopotentials were 
tested for transferability by comparing with all electron calculations for a range of 
typical atomic and ionic configurations. The pseudo-eigenvalues and total energies 
are equivalent to the all electron values to within a few me Vs. 

For Mn and O we constructed non-relativistic optimized pseudopotentials. The 
oxygen pseudopotentials were generated from a 2s 2 2p A reference configuration with 
core radii of 1.5 a.u. for both s and p orbitals. They were then optimized using 4 
and 3 basis functions with q c of 7.0 and 6.5 a.u. for s and p orbitals respectively. The 
oxygen pseudopotentials were used in earlier calculations for non-magnetic perovskite 
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oxides|T(J and gave accurate results. For Mn, a 4s°' 75 4p a25 3(f 5 reference configuration 



was used, with r c equal to 2.0, 2.15 and 2.0 a.u. for s, p and d orbitals. The 
optimized pseudopotentials for Mn were constructed using 4 basis functions and q c = 
7.5 a.u., and partial non-linear core corrections 

in 

were included. The core charge 



was approximated by a zeroth order spherical Bessel function within a radius of 0.737 
a.u., with the full core charge used outside this radius. The resulting optimized 
pseudopotentials are shown in Figures |] and |[ Again the pseudopotentials are 
transferable to within a few me Vs. 

All pseudopotentials were put into separable formpO] using two projectors for 



each angular momentum! 31 1. The first projector was taken as the atomic pseudo 
wavefunction, and the second projector as r 2 times the first projector, suitably or- 
thogonalized to the first |32||. For La, Bi and Mn, the I = component was chosen 



as the local potential. For oxygen the local potential was the I = 1 component. The 
absence of ghost states was confirmed using the ghost theorem of Gonze, Kackell and 
Schemer 

2.2 Technical Details 

The total energy and band structure calculations were performed on a Silicon Graph- 
ics Power Challenge L using the conjugate gradient program CASTEP 2.1 [p4|, [35| 
which we have extended to study spin-polarized systems. We used a plane wave cut 
off of 60 Ry, which corresponds to around 3500 plane waves in a cubic unit cell with 
lattice constant of around 4 A . A variable Gaussian broadening between 1 eV and 
0.002 eV was applied to the k-point sampling to speed convergence for metallic sys- 
tems. A 6x6x6 Monkhorst-Pack [pB{| grid was used for all calculations. This led to 10 
k-points in the irreducible Brillouin Zone for the high symmetry cubic structures, and 
a correspondingly higher number for the distorted structures with lower symmetry. 



The Perdew-Zunger parameterization ||37|| of the Ceperley- Alder exchange corre- 



lation potentially with the von Barth-Hedin interpolation formula |39l was used. 
Previous all-electron calculations for perovskite manganites||29|| found no appreciable 
difference between results obtained using the von Barth-Hedin and the more sophis- 
ticated Vosko- Wilk-Nusair |j40|| interpolations. 

For density of states (DOS) calculations we interpolated the calculated eigenvalues 
to a grid of ~ 350,000 k-points in the irreducible simple cubic Brillouin Zone using 
the interpolation scheme of Monkhorst and Pack [|36|. We then applied the Gilat- 
Raubenheimer method [f4~T| to this fine mesh. Finally, for the band structure plots, and 
for use in the tight-binding analysis, symmetry labels along high-symmetry lines were 
assigned using projection operators for the corresponding irreducible representations. 
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2.3 Comparison with LAPW calculations 



Before presenting our new results for BiMnOs, we first show that the results obtained 
for perovskite manganites using this PWPP implementation are in excellent agree- 
ment with the results of previously published LAPW calculations. Figure ^| shows 
the band structures for up- and down- spin electrons in cubic ferromagnetic CaMn03, 
calculated in this work using the PWPP method. We used a 10 electron Ca pseu- 
dopotential with a 3s 2 3p 5 3d 1 Ca 2+ reference configuration and core radii of 1.29, 1.29 
and 1.45 a.u. for the s, p and d potentials. Four basis functions were included in the 
optimization [Q, with a cutoff wavevector of 7.0 a.u. The band structures are indis- 
tinguishable from those shown in Figure 7 of Ref. [29], which were calculated using 
the LAPW method. In addition, we find that our PWPP calculated energy difference 
between the ferro- and paramagnetic phases of CaMnOs is 876 meV, within 2% of 
the LAPW value of 860 meV. The inclusion of non-linear core corrections is essential. 
The pseudopotential band structure calculated without non-linear core corrections is 
markedly different from the LAPW band structure. In particular the energy of the 
minority spin Mn 3d bands is too high, a pseudo-gap opens up in the majority spin 
O 2p bands, and the lower O 2p bands have unphysical flat regions. In addition 
the total energy difference between the ferromagnetic and paramagnetic phases is 
overestimated to be 3.2 eV. 



3 Calculated Electronic Properties of BiMnOs and 
LaMn0 3 

In this section we investigate the origin of the differences between BiMnOs an d the 
rare earth manganites by comparing the calculated electronic properties of BiMnOs 
with those of LaMnOs. We begin by calculating the electronic structure for the high 
symmetry cubic phases, without including magnetic effects, then lower the magnetic 
symmetry to the ferromagnetic phase. Finally we introduce structural distortions in 
both paramagnetic and ferromagnetic calculations. This ability to isolate structural 
and magnetic distortions is unique to computational studies, and allows identifica- 
tion of the essential microscopic interactions which cause the observed macroscopic 
behavior. All calculations are carried out at the experimental lattice constant of 3.95 

A. 
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3.1 Cubic Paramagnetic Structures 

First we present results for BiMnOs and LaMnOs in the highest possible symmetry 
state - that is the cubic structure, without spin polarization (we call this the param- 
agnetic (PM) phase). Although this phase is experimentally inaccessible, it provides a 
useful reference for understanding the ferromagnetic (FM) structures to be discussed 
later in this paper. 

Figure [7| shows the calculated densities of states for cubic paramagnetic LaMnOs 
and BiMnC>3. The plotted energy range is from -12 eV to 4 eV, and the lower lying 
semi-core states have been omitted for clarity. The Fermi level is set to zero in both 
cases. The broad series of bands between -2 and -7 eV in both materials arises from 
the oxygen 2p orbitals. Above the oxygen 2p bands, and separated from them by an 
energy gap, are the Mn 3d bands. The Mn 3d bands are divided into two sub-bands 
- the lower energy t 2g bands, and the higher energy e g bands - as a result of crystal 
field splitting by the octahedral oxygen anions. One striking difference between the 
two DOS plots is the presence of a band between -10 and -12 eV in the BiMnOs band 
structure which does not exist in the LaMnOs case. This band corresponds to the 
high lying occupied Bi 6s electrons. In addition, the high energy La 5d electrons have 
a very different form than the Bi 6p electrons, which occupy a similar energy range. 

The Fermi level lies near the top of the Mn 3d t 2g bands and is in a region of 
high density of states. The large DOS at the Fermi level suggests that the cubic PM 
structure is unstable, and that a lower energy structure could be achieved by allowing 
spin-polarization and/or structural distortion. 

Figure § shows the corresponding band structures along the high symmetry axes 
of the simple cubic Brillouin Zone. Again the broad O 2p bands between -2 and -7 eV 
can be seen clearly in both materials, with the Mn 3d bands above them, separated 
by an energy gap. The highlighted lines in the band structure plots accentuate the 
upper and lower Mn 3d bands, which have a similar form to each other and to those 
of CaMnOs (Figure |(|). This indicates a 'universality' in the manganite structure 
throughout the perovskite manganite series, which is independent of the identity of 
the large cation. 

Again we observe the Bi 6s band, which can be seen to have considerable dispersion 
suggesting that it is involved in covalent bonding. In addition, the Mn 3d bands in 
BiMnOs overlap with the partially occupied Bi 6p orbitals, whereas in LaMnOs the 
next highest bands are the unoccupied La 5c? bands. 

The differences between the two band structures show up clearly in the bands 
along the r — > X line. This region of the band structure is shown in Figure |9| with 
the symmetry labels added, and the energy scale extended to include the lower energy 
O 2s and La 5p bands. In the following discussion we focus on the bands with Ai 
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symmetry as capturing the essential features. The Ai bands are plotted with solid 
lines, and all other bands with dashed lines. 

In LaMnOs, the oxygen 2p Ai bands decrease monotonically in energy from T to 
X, and the Mn 3d A% band increases monotonically from T to X. Figure 10 shows 
a contour plot of the charge density in the group of bands above the O 2p manifold 
near the Fermi level at the T point. The spacing of the outer three lines is | that of 
the innermost contours. The charge density has been projected on the (100) plane 
by integrating through the whole unit cell perpendicular to (100). The Mn ion is 
located at the center of the plot, the La ions at the corners and the O ions in the 
middle of the cell edges and at the center. It is clear that the Fermi surface consists 
largely of Mn 3d electrons, with contributions from other atoms below the resolution 
of the plot. The charge density at the X point is very similar, except that there is a 
small oxygen contribution, indicating that the amount of Mn 3d - O 2p hybridization 
increases along the T to X line. 

In BiMnC>3 the behavior is quite different. First, the X\ symmetry of the Bi 6s 
band at the X-point causes the X\ O 2p band to be 'pushed up' in energy, resulting 
in a different ordering of the O 2p bands at the X point. The Mn 3d A x band crosses 
the very dispersive Bi 6p Ai band, and the latter moves below the Fermi level near 
the X point. Figures [11] and [12] show contour plots at the F and X points respectively, 
of the charge density in the five bands near the Fermi level above the O 2p manifold. 
(These are the bands which we found to be entirely Mn 3c? in the case of LaMnOs.) 
Again the charge density is projected onto the (100) plane, and the same scale is used 
for the contours as in Fig |10|. At the T point the charge distribution is similar to that 
of LaMnCV However at X there is a significant amount of charge density around the 
Bi atoms. Further analysis reveals that the Bi component is in the X' A band, which 
crossed the highest Mn 3d band. The lack of anti-crossing between the Bi 6p and Mn 
3d Ai bands is due to the zero matrix element between the Bi 6p and Mn 3d orbitals. 

To quantify the differences between cubic PM BiMnC>3 and LaMnOs, we per- 
formed tight-binding analyses of the T to X regions of the respective band structures. 
Tight-binding parameters were obtained by non-linear-least-squares fitting^] to the 
calculated ab initio energies at the high symmetry T and X points, and at 19 points 
along the A axis. First we performed analyses with only oxygen 2s and 2p and Mn 
3d orbitals included in the basis set. The tight-binding parameters thus obtained 
are given in Table |I[ The band energies calculated using these parameters have root 
mean square (RMS) deviations from the ab initio energies of 0.20 for LaMnOs an d 
0.25 for BiMnOs. The resulting band structures for the bands of Ai symmetry are 
compared with the ab initio values in Figure [13|. The limited basis set reproduces the 
LaMnOs bands well, consistent with an early proposal by Goodenough[Q that the 



9 



magnetic properties of LaMnOs are determined by the Mn 3d - 2p hybridization 
only Note that the fit to the lower energy O 2s bands is the least good - these bands 
are very close in energy to the La 5p bands which have not been included in the 
fit. The behavior of the BiMnC>3 Ai bands is less well reproduced, confirming that 
additional orbital overlaps are essential in producing the observed band structure. 

We then repeated the fitting procedure for BiMnOs, with Bi 6s and 6p orbitals 
added to the basis. Transfer integrals between nearest neighbor Bi 6s - O 2p, Bi Qp 
- O 2p, Bi 6p - Mn 3d and Bi 6p - Bi Qp orbitals were allowed to be non-zero. This 
significantly improved the quality of the fit to the ah initio bands (see Figure |14|), and 
reduced the root mean square deviation to 0.12. The values of the new tight-binding 
parameters are given in Table 0. The largest transfer integrals involving Bi are the 
Bi 6s - O 2p and Bi 6p - O 2p a interactions, with the magnitude of the cr-bonded Bi 
6p - O 2p interaction being approximately 30 % larger than that of the Bi 6s - O 2p. 
Also large are the Bi Qp - Bi 6p a interactions, which cause the Bi 6p Ai band to be 
pushed down below the Fermi level. 



3.2 Comparison with Experiment 

The conclusion from our tight-binding analysis is consistent with recent measure- 
ments showing enhancement of charge ordering in Bi-doped CaMnOs j|5]]. Rao and 



co-workers observed that the charge-ordered state in Bi .3Ca .7MnO3 persists to a 
higher temperature than that in Lao.3Cao.7Mn03. They explained their observations 
by noting that the electronegativity of Bi enhances Bi-0 hybridization and in turn 
reduces the amount of Mn-0 hybridization. 

Additional evidence in support of this phenomenon can be extracted from Ref. 
| pE6| , in which a temperature-composition phase diagram for Bi^j-CazMnOs, showing 
the charge-ordered transition and the Neel temperature, is plotted. Comparison with 



similar data for Lai_ x Ca x Mn03 (in for example Ref. [fl7 |) confirms that the charge 
ordered phase persists to higher temperature in Bi-doped CaMnOs than in La-doped 
CaMn0 3 . The ability to tune the positions of phase boundaries by substitution of 
Bi for the rare earth ions might prove valuable in optimizing materials properties for 
specific device applications. 



3.3 Cubic Ferromagnetic Structures 

Next, we present the results of calculations in which the high symmetry cubic struc- 
ture is retained, but the electrons are allowed to spin polarize. We find that introduc- 
tion of spin polarization reduces the energy of both BiMnOs and LaMn03 by around 
1 eV per unit cell compared with the paramagnetic case. The most important obser- 
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vation of this section is that the differences between BiMnC>3 and LaMnC>3 which we 
observed in the paramagnetic calculations persist into the ferromagnetic phase, the 
PM to FM transition introducing the same kinds of changes in both materials. 

Figure |l5] shows the calculated densities of states for cubic ferromagnetic LaMnOs 
and BiMnC>3. The majority spins are represented by the solid line on the positive 
y-axis, and the minority spins by the dashed line on the negative y-axis. In both 
LaMnOa and BiMnC>3, the down-spin Mn 3d band is split off from the O 2p band, 
and has a similar form to the corresponding paramagnetic band. The up-spin Mn 3d 
hybridizes strongly with the O 2p and there is no band gap for the majority carriers. 
The up-spin DOS at the Fermi level in LaMn03 is still high, indicating that the cubic 
FM state has a high energy. This is consistent with the fact that the lowest energy 
spin polarization in structurally relaxed LaMn03 is antiferromagnetic[fP| . The DOS 
at the Fermi level in BiMnOs is somewhat lower suggesting that the FM phase should 
be more stable in BiMn03 than in LaMn03. For both compounds, the Fermi level cuts 
through the very bottom of the down-spin Mn 3d bands, and the conduction band is 
occupied almost entirely by up-spin electrons. Materials in which the electrons at the 
Fermi level are 100% spin-polarized are known as half-metallics, and are considered 
desirable for use in devices such as spin-transistors. 

Again the most obvious differences between the two electronic structures are the 
presence of the Bi 6s band between -10 and -12 eV, and the contrasting forms of the 
Bi 6p and La 5d bands. 

Figure [16] shows the up and down-spin band structures for BiMnOs and LaMnOa 
along the high symmetry axes of the simple cubic Brillouin Zone. There are many 
similarities between the FM and PM energy bands, and our earlier conclusions re- 
garding the origins of the differences between BiMn03 and LaMn03 continue to be 
valid. 

First we examine the up- and down-spin LaMn03 band structures, and com- 
pare with the PM LaMn03 band structures to determine the changes which spin 
polarization causes in a 'conventional' manganite. The states which correspond to 
non-magnetic atoms are unchanged from the paramagnetic state, and are identical for 
up- and down-spin electrons. For example, the dispersion of the lowest O 2p band is 
identical for up- and down-spin, and for the PM phase. Also the La 5d bands, which 
were above the Mn 3d bands in the PM state are unchanged in form and energy. 
The characteristic perovskite manganite Mn 3d pattern which we remarked on earlier 
persists in the FM phase, appearing around 2 eV higher for the down-spin electrons 
than for the up-spin because of the exchange splitting. The up-spin Mn 3d and O 
2p bands are strongly hybridized and there is no gap between them. However, the 
down-spin Mn 3d are split off from the O 2p bands by a larger gap than in the PM 
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case. As a result, the Mn 3d bands occupy the same energy region as the unoccupied 
La 5d bands. 

Next we compare the up- and down-spin BiMn0 3 band structures. As expected, 
the non-magnetic Bi 6s band and the lower O 2p bands are identical for up- and 
down- spin. In fact the overall form of the down-spin O 2p bands is almost identical 
to that of the corresponding PM BiMnC>3 bands. However it is different from that 
of the down spin LaMnC>3 O 2p bands - in particular the structure at the T point 
is quite different, a consequence of the unusual behavior of the Ai band which we 
observed earlier in the PM bands. Both up- and down-spin Mn 3d bands show the 
characteristic perovskite manganite features (most noticeable in this case along the M- 
R-X directions) which we noted earlier in our discussions of the PM band structures, 
and of CaMnC>3. The form of the down-spin Mn 3d/Bi 6p bands is similar to that 
of the PM phase, although the energy separation from the O 2p bands is larger. In 
contrast, the majority bands show strong hybridization between the Mn 3d and O 2p 
electrons. The lowest Mn 3d bands at the T point (r' 25 ) move below the upper O 2p 
bands. At the R point, the Bi 6p bands are now above the Mn 3d bands, whereas in 
the down-spin and paramagnetic band structures they were below. 

The contrast between the BiMnC>3 and LaMnC>3 minority conduction electrons is 
striking. Figure is a contour plot of the small amount of down spin charge density 
in the region above, and separated by a gap from, the O 2p bands. The charge density 
through the entire unit cell has been projected onto the (100) plane. The electrons 
are localized entirely in the Mn 3d orbitals with no hybridization with the O atoms. 
This is a result of the large energy separation between the Mn and O bands, which 
in turn is a consequence of the minority Mn 3d bands being pushed to higher energy 
by exchange forces. There is no contribution from the La atoms, because the La 5d 
electrons never intersect the Fermi level. Figure [18] shows the charge density of the 
down-spin conduction band electrons in BiMnC>3. Rather than summing over the 
unit cell, two different slices through the unit cell are shown. The upper plot is in the 
Mn-0 (100) plane, with the Mn ion at the center and O ions at the edge mid-points. 
There is minority spin conduction band charge density localized on both the Mn and 
O ions. The lower plot is in the Bi-0 (100) plane, with Bi ions at the corners and an 
O ion at the center. We see that the Bi contribution to the charge density is both 
large and highly directional. The fact that the conduction electrons partly occupy 
p-type atomic orbitals should produce quite different transport characteristics than 
those observed in conventional rare earth manganites, where the conduction bands 
are entirely Mn 3d-type. 



12 



3.4 Ferroelectric distortions 



Finally we investigate the origin of the proposed ferroelectricity in BiMnOs. In keep- 
ing with the philosophy of this paper, we study the lattice distortions of the high 
symmetry cubic phases, and compare our results for BiMnOs with those for LaMnOs. 
We restrict our discussion to zone center phonons. 

The perovskite manganites have five atoms per unit cell, which results in 15 
phonon branches, 3 acoustical and 12 optical. At the T point all phonons are three- 
fold degenerate, so there is one acoustical phonon frequency (which is zero), and four 
optical frequencies. We are interested in the optical phonons which have negative 
eigenvalues, indicating lattice instabilities. 

The force constant matrices for LaMn03 and BiMnOs were determined by cal- 
culating the Hellmann-Feynman forces resulting from the displacement of each atom 
in turn 0.1 A along the z direction of the unit cell. The forces exerted on the Mn 
ions by the other ions were determined using the acoustic sum rule. We calculated 
the Mn-Mn force by applying the acoustic sum rule to both the columns and the 
rows of the resulting matrix. The two values differed by less than around 0.001 
Our calculated force constant matrices for paramagnetic LaMnOs and BiMnOs are 
tabulated in Table §. 

The paramagnetic cubic phases of both LaMn03 and BiMnOs have two unstable 
phonon modes. The frequencies and eigenvectors of the unstable modes are given in 
Table |j. The mode involving motion of the equatorial oxygens, shown in Figure |19] (b), 
does not correspond to a ferroelectric distortion and will not be discussed further. The 
displacement pattern for the second unstable mode is shown in Figure [B| (a). Here 
the large (Bi or La) cations are moving in opposition to the oxygen cage, resulting in 
a ferroelectric displacement. The imaginary frequency of this mode is twice as large 
in BiMnOs than in LaMn03, indicating a stronger instability in the Bi compound. It 
is interesting to note that the Mn is moving in the same direction as the oxygen ions. 
This is opposite to the behavior of the Ti ion in BaTi03 and PbTiOs, but similar to 
the behavior of Zr in PbZr03 [|16] . 



The imaginary frequency phonons in the ferromagnetic phases are less unstable 
than the corresponding imaginary frequency paramagnetic phonons. In fact the ferro- 
electric mode in LaMn03 is now only slightly unstable, at 21. li crrT 1 , and the other 
formerly unstable mode now has a positive frequency. However the frequency of 
the ferroelectric mode in ferromagnetic BiMnOs remains strongly unstable, at 82.30i 
cm -1 . 

It is clear from our analysis, that the presence of ferroelectricity in BiMnOs, and 
absence of ferroelectricity in LaMnOs, can be explained by the different zone cen- 
ter lattice instabilities. Although a definitive prediction requires calculation of the 
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phonon dispersion throughout the entire Brillouin Zone, it is likely that the weakly 
unstable zone center phonon in LaMn03 will be overshadowed by a stronger insta- 
bility elsewhere in the Brillouin zone, reproducing theoretically the experimentally 
observed Jahn- Teller distortion. Similarly, the very unstable ferroelectric mode at 
the zone center in BiMnC>3 is likely to dominate over possible unstable modes at 
other frequencies, confirming theoretically the suggested existence of ferroelectricity 
in BiMnC>3. We are in the process of developing spin-polarized density functional 
perturbation theory code|l7| which will allow us to investigate efficiently the nature 
of the phonons throughout the entire Brillouin zone. Our results will be the subject 
of a future publication. 

An alternative suggestion for the origin of the ferroelectricity in BiMnC>3 is that 
it originates from an electronic phase transition, in contrast to the conventional dis- 
placive ferroelectricity caused by a lattice distortion Such electronic ferroelec- 
tricity has been demonstrated theoretically {|9| in the insulating phase of the Falicov- 
Kimball model, and there is some experimental data which is consistent with the 
occurrence of electronic ferroelectricity |49|]. The correlations associated with the elec- 
tronic instability are included via a BCS-like wavefunction and thus the resulting 
energy is not expected to be well described within the local density approximation. 
In similar situations, however, LDA results have been successfully used to compute 
values of model parameters for input to many-body investigations |50|j . This is clearly 
a subject which merits further experimental and theoretical study. 



4 Summary 

Our results indicate that covalent bonding between the bismuth cations and oxygen 
anions in BiMn03 introduces additional orbital interactions compared with the rare 
earth manganites, in which the rare earth - oxygen interaction is essentially purely 
ionic. These additional orbital interactions in turn stabilize different magnetic and 
structural phases. Our results are consistent with the limited experimental data on 
doped BiMn0 3 -based compounds, and suggest that a modern experimental study of 
pure BiMn03 might reveal novel and potentially useful phenomena. 
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BiMn0 3 


LaMn0 3 


Eo2s 


-17.680595 


-17.836376 


Eo2p 


-3.444399 


-4.514022 


EMnZd 


-1.202430 


-1.191285 


V02s-02s 


-0.252108 


-0.243683 


V{02p-02p)a 


0.735211 


0.620568 


V(02p-02p)ir 


-0.132201 


-0.063466 


[V(02p-02p)a]2 


-0.211935 


0.182635 


W{02p-02p)-!r\2 


-0.044754 


0.082951 


Vo2s-Mn3d 


-1.720688 


-1.735814 


V(02p-Mn3d)a 


-1.964201 


-1.838490 


V(02p-Mn3d)iv 


1.036316 


0.878961 


V(Mn3d-Mn3d)S 


-0.003391 


0.066346 



Table 1: Tight-binding parameters for BiMn0 3 and LaMn0 3 obtained by non-linear- 
least-squares fitting to the ab initio eigenvalues along T to X. E indicates an orbital 
energy, and V an inter-atomic transfer integral. All transfer integrals are between 
nearest neighbors, except those with the subscript '2' which are between next nearest 
neighbors. Only the parameters listed in the table were allowed to be non-zero in the 
fitting procedure. 
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BiMn0 3 




-10.310688 




0.202695 


Eo2s 


-17.717958 


Eo2p 


-3.725773 


EMnZd 


-1.167871 


VsiQp— Bi6pa 


0.848443 


VsiQp— BiGpir 


0.166061 


VBi6s-02p 


-0.812502 


V(Bi6p-02p)a 


-1.061660 


V(Bi6p-02p)TT 


-0.145201 


V(Bi6p-Mn3d)a 


0.170850 


V(Bi6p-Mn3d)w 


0.078144 


Vo2s-02s 


-0.237506 


V(02p-02p)a 


0.652059 


V(02p~02p)-rr 


-0.125450 


\V{02p-02p)o\2 


-0.010196 


\V{02p-02p)ir]2 


-0.002797 


Vo2s-Mn3d 


-1.645508 


V(02p-Mn3d)a 


-1.926533 


V(02p-Mn3d)w 


0.959768 


V(Mn3d-Mn3d)8 


0.029347 



Table 2: Tight-binding parameters for BiMnC>3 obtained by non-linear-least-squares 
fitting to the ab initio eigenvalues along T to X. An expanded tight-binding basis 
including the Bi 6s and 6p orbitals was used. The labeling scheme for the parameters 
is the same as in Table [T. Only the parameters shown were allowed to be non-zero. 
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BiMn0 3 LaMn0 3 





Bi 


Mn 


Oz 


Ox 


Oy 




La 


Mn 


Oz 


Ox 


Oy 


Bi 


-5.47 


-2.37 


-0.73 


4.28 


4.28 


La 


-0.05 


-2.64 


-0.71 


1.70 


1.70 


Mn 


-2.46 


11.12 


-8.66 


0.00 


0.00 


Mn 


-3.07 


15.10 


-11.57 


-0.21 


-0.21 


Oz 


-0.64 


-8.57 


13.03 


-1.91 


-1.91 


Oz 


-0.72 


-11.53 


17.04 


-2.40 


-2.40 


Ox 


4.28 


-0.09 


-1.82 


-2.19 


-0.18 


Ox 


1.92 


-0.46 


-2.38 


-0.01 


0.91 


Oy 


4.28 


-0.09 


-1.82 


-0.18 


-2.19 


Oy 


1.92 


-0.46 


-2.38 


0.91 


-0.01 



Table 3: Force constant matrices in cubic paramagnetic BiMnOa and LaMnOs. 



BiMn0 3 LaMn0 3 



v (cm r ) 


Bi 


Mn 


Oz 


Ox 


Oy 


v (cm 1 ) 


La 


Mn 


Oz 


Ox 


Oy 


72.39i 
98.20i 


0.0 
-0.43 


0.0 
0.09 


0.0 
0.16 


i 

0.62 


i 

0.62 


49.04i 
44.69? 


0.0 
-0.59 


0.0 
0.22 


0.0 
0.21 


i 

0.53 


i 

0.53 



Table 4: Eigenvectors and eigenvalues of the dynamical matrix which correspond to 
the unstable phonon modes in cubic paramagnetic BiMnOs and LaMnOs. 
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Figure 1: The perovskite structure. The small B cation (in black) is at the center of 
an octahedron of oxygen anions (in gray). The large A cations (white) occupy the 
unit cell corners. 
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Figure 2: Optimized ionic pseudopotentials for La. The core radius is 2.6 a.u. for 
all angular momenta. Four basis functions were included in the optimization, which 
used a cutoff wavevector of 7.5 a.u. 
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Figure 3: Optimized ionic pseudopotentials for Bi. The core radius is 3.0 a.u. for 
all angular momenta. Four basis functions were included in the optimization, which 
used a cutoff wavevector of 7.5 a.u. 
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Figure 4: Optimized ionic pseudo-potentials for Mn. The exchange-correlation poten- 
tial contribution from the core charge has been subtracted, in accordance with the 
non-linear core correction scheme of Louie et al. [B5[]. The core radii are 2.0, 2.15 and 
2.0 a.u. respectively for the s, p and d pseudopotentials. Four basis functions were 
included in the optimization, which used a cutoff wavevector of 7.5 a.u. 
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Figure 5: Optimized ionic pseudopotentials for oxygen. The core radius is 1.5 a.u. 
The s pseudopotential was optimized using 4 basis functions and a cutoff wave vector 
of 7.0 a.u., and the p pseudopotential used 3 basis functions and a cutoff wavevector 
of 6.5 a.u. 
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Figure 6: Up- and down-spin band structures for CaMnC>3 calculated using the plane 
wave pseudopotential method and the partial non-linear core correction method of 
Louie et al.[[2J|- The band structures are in good agreement with earlier LAPW 



calculations 29 
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Figure 7: Calculated densities of states for cubic paramagnetic LaMn0 3 and BiMn0 3 . 
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Figure 8: Calculated band structures for cubic paramagnetic LaMnOs and BiMnOs 
along the high symmetry axes of the Brillouin Zone. The highlighted lines in the band 
structure plots accentuate the upper and lower Mn 3d bands, which have a similar 
form to each other and to those of CaMnOs (Figure |). 
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Figure 9: Calculated band structures for cubic paramagnetic LaMn03 and BiMn03 
along the r — > X axis. The solid lines show the bands of Ai symmetry. 
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Figure 10: Contour plot of the charge density distribution in the Mn 3d bands at 
the T point of LaMn0 3 . A (100) projection is shown, with the Mn cation at the 
center of the square, the La ions at the corners, and the oxygen anions at the edge 
mid-points and the center. The charge density has been integrated through the unit 
cell perpendicular to the (100) plane. 
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Figure 11: Contour plot of the charge density distribution in the Mn 3d bands at 
the T point of BiMnC>3. A (100) projection is shown, with the Mn cation at the 
center of the square, the Bi ions at the corners, and the oxygen anions at the edge 
mid-points and the center. The charge density has been integrated through the unit 
cell perpendicular to the (100) plane. Note the similarity to the corresponding plot 
for LaMnOs (Figure 0). 
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Figure 12: Contour plot of the charge density distribution in the 'Mn 3cf bands at 
the X point of BiMn0 3 . A (100) projection is shown, with the Mn cation at the 
center of the square, the Bi ions at the corners, and the oxygen anions at the edge 
mid-points and the center. The charge density has been integrated through the unit 
cell perpendicular to the (100) plane. 
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Figure 13: Comparison of ab initio A\ bands with those obtained from a tight-binding 
fit using only Mn and O orbitals in the basis. The fit for LaMn03 is good indicating 
that only the Mn and O ions are significantly involved in covalent bonding. The 
fit for BiMnOs has a higher RMS deviation, and in particular misses the additional 
curvature around -3 eV near the X point. This shows that additional interactions are 
present. 
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Figure 14: Comparison of ab initio A 1 bands with those obtained from a tight-binding 
fit including Bi 6s and Qp orbitals in addition to Mn 3d and O 2s and 2p. The fit is 
improved over the fit using only Mn and O orbitals. 
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Ferromagnetic densities of states 
in cubic BiMn0 3 and LaMn0 3 
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Figure 15: Calculated densities of states for cubic ferromagnetic LaMn0 3 and 
BiMn0 3 . 
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Figure 16: Up- and down-spin band structures for cubic LaMnC>3 and BiMnC>3 along 
the high symmetry axes of the Brillouin Zone. 
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Figure 17: Minority conduction band charge density in cubic FM LaMnC>3. The 
charge density is projected on the (001) plane of the cubic unit cell. 
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Figure 18: Minority conduction band charge density in cubic FM BiMn0 3 . The 
upper plot is a (100) slice through the Mn-0 plane, and the lower plot is a (100) slice 
through the Bi-0 plane. 
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Figure 19: Eigenvectors of the two unstable T point phonon modes in BiMn0 3 . The 
Mn ion is at the center of the unit cell surrounded by an octahedra of oxygens, with 
the Bi ions at the unit cell corners. Each mode is three-fold degenerate, (a) is the 
mode which leads to a ferroelectric distortion; (b) is a non-ferroelectric mode. 
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